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FOREWORD 


This document is the final technical report summarizing work 
performed by the Microwave Division of Aeroj et- General Corporation, 

El Monte, California, during Phase II of the "Microwave Radiometer Data 
Study." Phase II of this study, sponsored by the NASA Manned Spacecraft 
Center, Houston, Texas, under Contract No. NAS 9-9363 is concerned 
with atmospheric and antenna pattern effects on measurements with the 
NASA P 3 A Multifrequency Microwave Radiometer System. This report 
summarizes equations and data flow developed to correct measurement 
data for these effects and includes results obtained from flight experiments 
designed to substantiate the equations. This work was performed between 
7 October 1969 and 15 May 1970. 

Personnel who have contributed to the data study and preparation 
of this report include A T Edgerton, A Stogryn, G Poe, S. Baker, 

G. T C half in, J. Kirwan, P. R Jordan and F Ruskey 
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Section 1 


INTRODUCTION AND SUMMARY 


This document is the final technical report summarizing 
work performed by Aerojet-General Corporation during Phase II of the 
Microwave Radiometer Data Study. Phase II of this study, sponsored 
by the NASA Manned Spacecraft Center under Contract No. NAS 9-9363, 
is concerned with the effects of atmospheric attenuation and emission, 
and the antenna patterns associated with the NASA P3A Multifrequency 
Microwave Radiometer System (MFMR). The object of this study is to 
develop equations and data flow for a computer program to correct 
measured radiometric data for these effects and to develop flight ex- 
periments and test procedures for verifying the equations. Data from 
this study were supplied to Manned Spacecraft Center’s Computation and 
Analysis Division for computer programming. 

This effort represents a portion of the over-all task of obtaining 
absolute brightness temperatures from raw MFMR data. 

Absolute radiometric brightness temperatures can be derived from 
MFMR data by correcting raw data for effects of 1) instrumentation 
calibration (baseline drift, gain, radiometer system losses including 
waveguide and antenna insertion losses, etc.), 2) antenna patterns, 

3) radome effects, and 4) atmospheric effects. The basic data reduction 
steps necessary for applying instrumentation and gross radome corrections 
to raw MFMR data have been developed by NASA/MSC and contractors 
The task of correcting measurement data for radome effects is the objective 
of a separate study. 

This report deals with techniques for removing atmospheric and 
antenna pattern effects from MFMR data Section 2 contains a detailed 
description of equations, data flow, and experimental methods for removing 
atmospheric effects from measured data Section 3 contains equations, 
data flow, experimental and data reduction procedures for correcting MFMR 
data for antenna pattern effects Section 4 is a discussion of verification 
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tests designed to confirm the data reduction equations and rationale. 

Section 5 contains conclusions and recommendations based on the overall 
Phase II study. 

In summary. Phase II of the Microwave Radiometer Data Study 
provides techniques for removing atmospheric contributions and antenna 
pattern effects from MFMR data The techniques are sufficiently general 
for use with other microwave radiometer systems, and computer programs 
are being developed by the Manned Spacecraft Center to apply these corrections 
to data acquired with the MFMR, Flight experiments were formulated and 
conducted to substantiate the data correction routines Data obtained 
during these experiments v/ere inconsistent. Unrealistic brightness 
temperatures were obtained using antenna and radome loss values furnished 
for use m this study. Significant improvement was achieved by using an 
empirically determined loss factor However, quantitative results could 
not be obtained due to the sensitivity of measurement data to the empirically 
determined loss factor. After an adequate calibration of the MFMR is 
achieved, verification tests of the atmospheric and antenna pattern routines 
can be completed. The microwave radiometer data processing routines 
developed under Phase II of the data study will, for the first time, enable 
routine application of atmospheric and antenna pattern corrections. 

These refinements will facilitate use of MFMR data by the scientific 
community. 
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Section 2 

ATMOSPHERIC ATTENUATION AND EMISSION STUDY 


This section of the report documents equations and procedures 
developed to correct MFMR data for atmospheric attenuation and 
emission effects. 

2. 1 INTRODUCTORY COMMENTS 

The brightness temperature of a surface observed by a 
radiometer at some height z above the surface is affected by the 
atmosphere in several ways. First, the radiation leaving the surface 
contains a component which consists of reflected incoming radiation 
emitted by the ionosphere above the surface (plus a component due to 
the cosmic background). Secondly, on leaving the surface, the radiation 
is attenuated by the atmosphere in propagating to the radiometer. 

Finally, the atmosphere between the radiometer and the surface itself 
emits radiation which contributes to the received signal. These contri- 
butions to microwave radiometric measurements can range from a 
negligible amount for a clear, dry day at frequencies below 15 GHz, to very 
substantial values for humid, cloudy and/or rainy conditions at frequencies 
above 20 GHz. Thus, the atmospheric contributions to airborne 
microwave radiometric data must generally be removed to obtain 
the true brightness temperature of the terrain surface The task is to 
devise a method fox' removing atmospheric'contributions from radiometric 
measurements where the surface of interest and the intervening atmosphere 
are not necessarily uniform. 

MFMR missions are flown for a variety of purposes and con- 
sequently, requirements for atmospheric corrections vary from mission to 
mission. Flights are conducted at various altitudes and under a wide 
range of climatic conditions. To accommodate the range of experimental 
conditions, a flexible atmospheric effects correction method must be 
devised. Moreover, planning of MFMR missions must include consideration 
of data requirements and mission constraints for performing an adequate 
atmospheric correction. Uniform atmospheric conditions are generally 
desirable. 
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Atmospheric contributions to MFMR measurements can be 
determined, by (1) computation of atmospheric attenuation and 
emission using suitable mathematical expressions and known or estimated 
meteorological conditions, and/or (2) direct MFMR measurements of 
atmospheric radiation at selected heights above the terrain- Each 
approach possesses advantages and limitations. Both appear to be 
useful for MFMR data analysis. In the following discus sion both approaches 
are considered. 

To compute the effect of the atmosphere on measured brightness 
temperature, ideally one would like an air temperature and water vapor 
profile of the atmosphere together with a specification of the water 
droplet distribution within any clouds that may be present. Except for 
areas located near weather stations, such a detailed description of the 
atmosphere will not be available. Fortunately, the number of parameters 
that must be known to accurately predict atmospheric effects is small if 
the observations are made at sufficiently great height, (1 e. , above 20, 000 
feet). In fact, detailed numerical studies at a number of frequencies 
have shown that m the absence of ram, only the total water vapor content 
(the integral of the water vapor density, from sea level to the top of the 
atmosphere), the total cloud water droplet content, and the air temperature 
at cloud height will suffice to predict atmospheric effects to within I or 
2°K. The PSA aircraft possesses sensors for determining height 
above terrain, barometric pressure, relative humidity and air temperature, 
so that many of the needed parameters can be obtained from on-board 
sensors. Section 2. 2 contains a comprehensive discussion of numerical 
models for estimating atmospheric contributions to MFMR data. 

The direct measurement approach requires more flight time than 
that needed for numerical approximations of atmospheric contributions. 
However, better accuracy can be realized with direct measurements when 
nonumform weather, such as variable cloud conditions, is encountered. 

The number of MFMR sky measurements (flight altitudes) necessary to 
define atmospheric contributions to terrain observations is a function of 
prevailing weather conditions. Section 2, 3 contains a comprehensive 
discussion of the direct measurement approach. 
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2.2 


THEORY 


The brightness temperature of a surface observed by a 
radiometer at some height z above the surface is affected by the 
atmosphere in several ways. First, the radiation leaving the surface 
contains a component which consists of reflected incoming radiation 
emitted by the atmosphere above the surface (plus a component due to 
the cosmic background). Secondly, on leaving the surface, the radia- 
tion is attenuated by the atmosphere m propagating to the radiometer 
Finally, the atmosphere between the radiometer and the surface itself 
emits radiation which contributes to the received signal The result 
of the last two effects on the brightness temperature of the radiation 
leaving the surface can he expressed in a simple way if scattering effects 
in the atmosphere are ignored while the mathematical expression of 
the first effect can involve rather complex equations if the surface is 
treated in a precise manner.' 1 

Consider the situation illustrated in Figure 2-1 If the brightness 

temperature of the radiation with polarization p leaving the surface in 

the direction 6 is T (0, o), then the brightness temperature at the 

P 

height z is 

T p (6,a) = L(6,z) T p (G,o) + T atm (6, Z ) (2-1) 


where the loss factor L is given by 


li(8, z) = exp 


-sec 8 


J* x(u)du 


(2-2) 


and the temperature of the radiation emitted by the atmosphere between 
the ground and the radiometer is 


T atm ( 2 ) = sec 8 J T a - r (z 7 ) ^(z 7 )exp -sec 8jx(u)du 


dz ' 


(2-3) 


For the frequencies of interest, scattering in the atmosphere can be 
ignored except for heavy ram storms 
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FIGURE 2-1 


GEOMETRY FOR ATMOSPHERIC CORRECTIONS 
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In (2-2) and (2-3), x(z) is the absorption coefficient per unit length m 
the atmosphere at the height z while T . (z) is the air temperature 
at the height z Note that X is a function of the radiometric frequency 
although, in order to simplify the notation, the frequency dependence is 
not indicated in the equations 

As mentioned above, the expression for Tp(0, o) m terms 
of the properties of the surface and the atmosphere may be very complex 
Assuming that the incident sky radiation is unpolarized and that the 
ground temperature does not vary significantly over the depth from 
which most of the radiation is emitted, 

T p( e, o) = V 0)T g + ?fJ[Vph + V] hky < e '> dfi ' <^) 

where e (9) is the emissivity of the upper hemisphere ground in 

ir 

the direction 9 , the thermometric temperature of the ground 
an< ^ ^ph a-nd y a re respectively coefficients for scattering with 
polarization p from the ground when the incident radiation has 
polarization h (horizontal) and v (vertical). Both y^ and y 

depend on the incidence angles 0 ' , and scattering diiection 0, <? 

The quantity T g ^ ( 6 ') is the sky temperature which is given by 

T sky* 9 ^ = sec 8 J T air ( z O^z^exp [-sec 6 ' Jx(u)du ]dz ' + T CB L(8'eo) (2-5) 

O 1 O j 

In (2-5), T air > x , and L have the same meaning as above while 

2 7 K) is the temperature of the cosmic background radiation 
The convention that z =<o corresponds to the top of the atmosphere is 
used in (2-5) although in practice a finite number of the order of 20 km 
can be used to give adequate accuracy m numerical computations It is 
readily apparent, even apart from the problem of specifying the scattering 
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coefficients of the surface, that the use of (2-5) in (2-4) can lead to 
computations on a rather large scale. Therefore, the assumption is 
often made that the surface under investigation is perfectly flat and 
uniform (specular surface). In this case, the scattering coefficients 
contain delta functions of angle and (2-4) reduces to 

Tp (0,o)= V 0 )T g + [l - y e >Kky(9) 


For many applications, the set of equations (2-l), (2-2), (2-3), (2-5), 
and (2-6) provide a useful means of accounting for atmospheric effects 
even if the surface under observation is not strictly specular 


2 2 1 Basic Integrals 

The preceding discussion has indicated that atmospheric 

effects can be expressed in terms of the three basic quantities L(9, z), 

Tat m (®» z ) anc * given by (2-2), (2-3) and (2-5). For realistic 

models of the atmosphere, it is clear that the integrals which occur m 

the definitions of these quantities cannot be evaluated analytically, 

and hence, a numerical quadrature must be performed. The most 

efficient quadrature scheme is obtained by noting that, except for the 

factor T . (z '), the integrand in both (2-3) and (2-5) are exact 
air 

differentials. Thus, if the integral over z ' in these equations is written 

as a sum of integrals, each over a sufficiently small interval m which 

T a j r (z ') may be considered a constant, the integral over z / may be 

performed analytically. For example, if z corresponds to the top of 

the atmosphere, and the interval (O, z ) is subdivided by the points 

mELX 

0=z.<z- ) <z 2 .. <z = z in such a way that T (z) is effectively 

1 2 n max 3 air 3 

a constant T (z.) m the subinterval z. < z < z , , , then 
air' V i i+l’ 


z N-l Z i+1 z 

(0) = T^^exp |-sec 0j*x(u)duj+ X T a - r (z.) f sec 0 x(z Oexpj-sec sj*x(u)du]dz 

o i=l *is. o ' 



r 2 i^ 1 - 

[ r r z ] 

= T CB exp 

-sec 0J X(u)du +_£ T air (z.) 
o J i= 1 

1 -exp | -sec 0) X(u)du 
1 *• o 


i 

z =z 


i+l 


/ 

z =z 
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(2-7) 



N-l 


z 

+ Z T (z ) 
air x' 

exp 

{* 1 

-sec 0 J x(u)du 

1=1 


o 


-exp 


[- 


sec 




If the height of the radiometer, z , is chosen to coincide with one of 
the above division points {say /), then a similar calculation shows 
that - 


N-l 

T atm< 8 ’ Z N'^-£ T air< Z i> 
1=1 


- z n' 


exp [- sec •P x (u)duj -exp [-sec 6 J' x( r a)du|J (2- 

rr ^ 


8 ) 


l+l 


Comparing (2-2), (2-7), and (2-8), it is seen that if the integral 


X z 

1 H(u)du 


(2-9) 


is calculated for l = 1, 2, . , N, then the three basic quantities can be 

expressed as 


3->{ e, z N /) = exp 


-sec 0f(z N /)] 


( 2 - 10 ) 


T 


atm (0J 

T sky 


N-l 


Z N ') 

= z 
1=1 

T air (z i ) j ex P 

(8) 

- 

T^^exp^-sec 


[-sec 6 (f(z N /)-f(z i+1 ) 

-exp [-sec 9 [j 

:(Z N *-«*!>)] 

N-l 


(2-11) 

ef(zj] +r t (z ) 
' N 7 J . , air v v 

i=l 

| exp [-sec Sf(z^) 


- exp [-sec 0f(z i+1 )jJ 

(2-12) 
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To evaluate these expressions, all that remains is to specify the 
dependence of the absorption coefficient x on atmospheric conditions 
and radiometric frequency. 

2. Z. Z The Absorption Coefficient 

The absorption coefficient X (s) may be written as a sum 

X (z) = x,_ (z) + x_ T ~(z) + x Jz) + x (z) (2-13) 

' 1 H 2 0' cloud' ' ram' ' ' ' 

where X_ , X^ q, X cloud’ an ^ X rain are res P ect ively the absorption 

coefficiemrs of oxygen, water vapor, condensed water droplets (liquid 

or ice) m clouds, and ram. The contribution of X . to X will not 

be considered further here because the occurrence of rain is somewhat 

less common than the occurrence of cloud cover and because X depends 

ram 

strongly on the rain drop size distribution function as well as rainfall 

rate, and hence, its computation depends on a detailed knowledge of the 

characteristics of a given rain storm The other three terms in (2-13) 

will be discussed below. The coefficients x „ and x, T 

O,, H^O are 

usually expressed in terms of the attenuation y i- n db/km_ The 

relation between y and X is 

X = .2302581 Y km" 1 (2-14) 


2. 2. 2. 1 Oxygen Attenuation 

The theory of absorption by oxygen was first given by Van 
Vleck [2-1] . Substituting numerical values for the coefficients 
and assuming an atmospheric content of 20. 9% oxygen by volume, 
Van Vleck r s formula for absorption by oxygen m air becomes 


= 2 6742 -4^- v2 ^ 
u 2 T . odd K 

air 


K(2K+3) 0 ^ (K+l)(2K-l) c 

“K+I S K + K S K- + 

T 


2(K‘ 2 +K+1)(2K+1) „ 
KT-KTT) b K 


exp 


-2 06844K(K+1)/T 


air 


db/km 


(2-15) 
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where 


P 


air 


T . 
air 


v 


= air pressure in mm Hg 
= air temperature m °K 
= i adiometric frequency in GHz 
AV , AV 


K± 


< V ~ V K ±) 2 + < AV ) 2 (" J+V K =)^ + (AV ) r 


Av 


Ko v ? ‘ )• (Av) ? ' 


AV = line width constant m GHz 


= re sonant line frequencies in GHz 
The first 46 resonant frequencies are listed in Table 2-1. 

Experiments have shown that, for all lines except the 1- line, 
the line width constant is given by [2-2] 


AV = 1 9xl0 3 P . I 760 
air 


, , . 290 V 9 GHz 

V 760-hP . )\ T / 

\ air J \ x air / 


The 1- line width constant has been determined as 1. 90 MHz /mm Hg 
at 300°K [^2-3j and is proportional to the air pressure Assuming that 
the width varies as the reciprocal of the air temperature, this leads to 


Av i_ = *- 444 


300 


air , 


GHz 


(2-16) 


(2-17) 
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Tattle 2-1 


MICROWAVE LINE FREQUENCIES (IN GHz) OF 0 2 


v l+ 

rr 

56 264752 

V l- 

= 

118 7507 

V 3+ 

- 

58. 446577 

V 3- 

=r 

62. 486256 

V 5+ 

= 

59* 590979 

V 


60 306055 

V 7+ 

- 

60. 434779 

V 

= 

59, 164205 

V 9+ 

- 

61. 150567 

V 

= 

58. 323885 

V ll+ 


61. 800163 

v n- 

= 

57. 612505 

V 13+ 

= 

62. 411224 

V 13- 

= 

56. 968245 

V 15+ 

= 

62 997980 

V 15- 

= 

56. 363448 

V 17 + 

- 

63 568535 

V 17- 

- 

55 783884 

V 19+ 

- 

64 1276 

V 19- 

= 

55 221449 

V 21+ 

= 

64. 6782 

V 21- 

= 

54. 6716 

V 2 3+ 

= 

65 223961 

V 23- 

- 

54. 1309 

V 25+ 

- 

65. 7626 

V 25- 

- 

53. 595951 

V 27+ 

" 

66. 2978 

V 27- 

= 

53. 0695 

V 29+ 

= 

66. 8313 

V 29- 

= 

52. 5458 

V 31+ 

= 

67. 3627 

V 31- 

= 

52. 0259 

V 33+ 

= 

67. 8923 

V 33- 


51. 5091 

V 35+ 

= 

68. 4205 

V 35- 


50. 9949 

V 37+ 

- 

6 8 9478 

V 37- 


50. 4830 

V 39+ 

= 

69 . 4741 

V 39- 

= 

49. 9730 

V 41+ 

- 

70. 0000 

V 41- 

= 

4-9. 4648 

V 43+ 

= 

70. 5249 

V 43- 


48 9582 

V 45+ 

- 

71. 0497 

V 45- 

= 

48. 4530 
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2 2. 2. 2 Water Vapor Attenuation 

Water vapor has absorption lines located at 22 33, 183 6 
and 324 GHz, as well as a great number of lines in the submillimeter 
range The attenuation due to water vapor can be written as 


Y H ? 0 ~ Y H 0(22) + Y H ? 0(183) + y H ? 0(324) + Y H,0(>324) 


(2-18) 


All numerical coefficients which are given below are based upon ex- 
perimental measurements. 

For the 22 GHz line [2-4] 

2 


% 0(22)= 4 - 23xI ° 2 A /S P H z O 

A air 


where 


Av 


Av 


(V-22. 33) 2 + (AV) 2 (V+22 33 2 +(Av) 2 


|db/km 


(2-19) 


AV 




( 2 - 20 ) 


Here and in the following 

3 

P H o = water vapor density in gm/m 
2 

T = air temperature in °K 
air 

P^. = air pressure m mm Hg 

v s= radiometric frequency m GHz 

For the 183 GHz line [2-5] 


y K, 0(183) = *• 361x10 ~ S/2 P H,0 


T air 


where 


v 


AV 


(v-183. 6r+(Av) tf (V+183. 6) <s +(Av) i 


Av=5 - 4 (w-) (|~) 


[db/km 
( 2 - 21 ) 

( 2 - 22 ) 
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There is less reliable information on the 324 GHz line. We 


will write 


0(324) - 4 22x1 ° 2 P H„0 


air 


[ 


Av 


+ 


Av 


2 (v- 324) + (Avr (V+324r+(Av)' 


db/km (2-23) 


and assume that 


* ■ 6 <T&-) (F- 

air, 


(2-24) 


Finally, for frequencies below 300 GHz, the effect of all 
higher frequency lines can be combined and written as 


2 

P . 

% O ( >22 4) = 1 510 — p^. ^ db/km 

2 T air 2 


(2-25) 


2. 2. 2. 3 Absorption by Condensed Water 

Clouds consist of minute droplets of water or ice If these 
droplets are assumed to be spherically shaped, the Mie theory of 
scattering is applicable. Since the paiticle sizes are so small compared 
to the wavelengths under consideration, the limiting form of the theory 
maybe used In the limit of small particle sizes, it can be shown that, 
after averaging over-particle sizes, the absorption depends only on the 
mass of the water or ice. The result is [2-6] 


X _ 6irx 10' 

cloud c 


v 



(2-26) 


where 

M = mass of droplets per unit volume 
p = density of water (or ice) 

K = dielectric constant of water (or ice) 
c = speed of light ( = 2. 997925xl0 8 m/sec) 
V = radiometric frequency m GHz 
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and Im means "imaginary part" of the bracketed quantity. 


Since the imaginary part of the dielectric constant of ice is 
small, it is found that ice particles contribute very little to the absorption 
coefficient unless they are near a melting condition in which case their 
surfaces may be wet. Consequently, attention will be confined to 
water droplets. Inserting numerical values for the constants, (2-26) 
becomes 


K . , = 6.2875342 x 10 

cloud. 


Mv Im 


[&*■] 


(2-27) 


where M is the density of the condensed water m gm/m and V the 
frequency in GHz. 

The dielectric constant of water, which appears m (2-27), 
may be expressed by a formula of the Debye type* 
e - e 

O 00 

K = r-m + e « 

1 -1 


(2-28) 


where e , , 


and the relaxation wavelength AX are functions of the 
temperature. A numbei of investigations of these parameters have 
been made. Hasted [2-7] gives 


e = 87.44 - . 4008 (T -273) 4 9. 398 x 10“ 4 (T -27 3) 2 4 1 41 x 10" 6 (T -273) 3 

■* Tr w w 

(2-29) 


w 


w 


where T is the water temperature in U K. This may be assumed equal 


w 


to the air temperature T . at the height of the cloud There is less 
information about the temperature dependence of c ro . The current 
evidence suggests that it is independent of temperature. We will take 
When v is expressed in GHz, is given by 


£ oo = 5.5 


10 

~n 


-4 


[-. 555(T w -263)(T w -253)(T w ~ 243) + 1. 19 (T w -273) (T w -253) (T w -243) 


-.87 (T -273) (T -263) (T -243) + .227 (T -273) (T -263) (T -253)1 
' w w w w w w J 


(2-30) 
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2. 3 DETERMINATION OF ATMOSPHERIC CORRECTIONS 


The factors determining the atmospheric contribution to the total 
brightness temperature of a target have been discussed m the previous 
section In determining the correction to be applied m a given situation, 
there are e ssentially two methods of approach. The first consists m 
assuming that atmospheric conditions are sufficiently well known (e. g. , 
from data supplied by a nearby Weathei Bureau station or by direct 
meteorological measurements from the aircraft) that an a priori 
calculation of the atmospheric effects can be made. This will be called 
Option A. The second approach (Option B) assumes that the necessary 
parameters for an a priori calculation are not known (or known only 
partially) and thus requires that some measurements be performed before 
atmospheric corrections can be made. Option B, of course, is a more 
general procedure than Option A, but has the concomittant disadvantage 
of requiring some flight time to be set aside specifically for atmospheric 
measurements. In addition, precise calibration of the radiometers is 
mandatory if the direct measurement method is to succeed. 

These two approaches are discussed m more detail below 

2. 3. 1 COMPUTATION OF ATMOSPHERIC ATTENUATION AND 

EMISSION - OPTION A 

If the air temperature, air pressure, atmospheric water vapor 
profile, and condensed water distribution m the atmosphere are known, 
the calculation of the three basic quantities L(8,z), T a j. m (6, z), and T g ^, 
which were discussed m the theory, proceeds m a straightforward 
manner. The computer program or subroutine will require the following 
input data 

A. Radiometric frequency or frequencies v (m GHz). 

B. Radiometer look angle (or angles) 0. 

C. Aircraft altitude z^ /(m kilometers), z^/ should 
coincide with one of the z mentioned m D. 

D. A description of atmospheric conditions m the form of a 
table This should consist of a sequence of altitudes 

z^ (in kilometers) l = 1, . . , N at which the air 
temperature (°K), air pressure (mm Hg), water vapor 
density (gm/m ), and cloud water droplet content 

O 

(gm/m 0 ) are specified. 
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The calculation proceeds m the following manner* 


A 


B 


C 

D„ 

E. 


F. 


At each z^, and frequency V of mtere st calculate Yq 

according to equations (2-15) - (2-17) Change y 
to Xq according to (2-14). 


At each z^, and frequency v of interest calculate 
Yjj q according to equations (2-18) - (2-25). Change 

z 

Y to x^ q according to (2-14) 

At each z^ and frequency of interest calculate * c i ou{ j 
according to (2-27) - (2-30). 

Calculate X according to (2-13). 

Evaluate the sequence of integrals (at each frequency v). 
z 

KzJ = £ \du i=l,„.„N 
by any standard numerical integration rule. The 
trapezoidal rule will yield adequate accuracy and 
provide an easy implemental scheme if z ~ z <,1/3 km. 
For each angle 0 and frequency v calculate L, T atm > 
and T g ^ according to (2-10) - (2-1 2). If the trapezoidal 
rule is used m E above, T ^{z ) may be taken simply as 


2 


T (Zil) + 
air l+l 


T (z.) 
air r 


2. 3. 2 DIRECT MEASUREMENT OF ATMOSPHERIC EFFECTS 
OPTION B 


If sufficient a priori information concerning atmospheric 
conditions is not available, then specific measurements for this purpose 
must he made. The basic set of quantities that need to be estimated for 
atmospheric corrections is Li(0, z), T a ^ m (0, z) and T g j as m Option A, 
It will be assumed that the angle 0 is fixed m the following discussion 
(recall, also, that these quantities are frequency dependent, although the 
frequency dependence has been suppressed m the notation). 


The required measurements consist of the following. 

Starting at a very low altitude (i. e. , sufficiently low so that the atmosphere 

between the aircraft and ground does not contribute significantly to L, 

T , , or T , - typically this may be about 300 to 500 feet) and with 

3-trn sicy 
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the radiometers pointed skyward, a measurement of the brightness 
temperature is made at each of the frequencies of interest. A measurement 
of the air temperature is also necessary. The aircraft then increases 
its altitude by an increment Az and the measurements are repeated. This 
process is continued until the aircraft has reached an altitude at which 
the measured brightness temperature does not differ significantly 
from the cosmic background temperature or, if this is not feasible, to 
an altitude above which little water exists m the air {the contribution of 
the atmosphere above this altitude to the three basic quantities may be 
approximated with sufficient accuracy by the use of an a priori model). 


The data which is so obtained may be analyzed m the following 

manner. At the lowest altitude, z^, the brightness temperature is 

appi oximately the quantity T g ^ . The loss factor L(9, 

equal to 1 and z^) is approximately equal to zero (see equations 

(2-2) and (2-3) in Section 2.2). These estimates may be used to determine 

L(9, z) and T a ^. m {6, z) inductively. Thus, suppose that at some altitude 

z^, L/{6, z^) and T^_^(Q, z^) are known To determine L(0, z^+Az) and 

T , (9, z +Az) we may write equation (2-2) as 

Etm l 



rz + Az 

T (0,z+Az) = sec 0 1 T (z) K(z )exp 

atm' * i ' J air' ' ' ' ^ 


•z f Az 


1 


-sec 0 /; 1 K(u)du 


dz 


= sec 0 1 


10 1 


z + Az 
i 




( 2 - 


+ Az 


= e 


-sec y 


f 1 K(u)du ^ 


atm 


( 0 , 


z i> 


n Z + AZ 

z + Az 

+ sec 9j 1 T <«') K(»0«P 

-sec 8 J Kdu 

z 1 | 

z 


But, if the altitude increment Az is not too large, then T (z ') may, 
under most circumstances, be well approximated by T 


<*,) = 7 [ 

U1A JL C-t [_ 

T (z^ + Az)J will be a known quantity. Hence, we see that 


T (z ) + 
air' r 
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L(0, z + Az) 


T Mm (e > a i + Az > = TTiTin ^{8,* ) + ?,,>,) 1 - T 7 B V T • (2-32) 


atm 1 5 1 ' 1 ‘air' V 


L(0,z -I- Az) 


At the same time, the difference m the measured brightness temperatures 
at altitude s z. and z^ + Az is given by 

T b ( e » z i) " T b (9, z^Az) = s ec 6 j T air^ Z ^ K ^ ^ exp “ sec 6 / Kdu dz * 

z z 

i L i 


+ -mv) T CB - T b< e '\ +Ai > 


■ z + Az m 

f 1 . f + 

J + j + 

z z + Az 1 

L i l J 


= sec 8 


I#? T CB -T b (6, 2i+ Az). 


r z,+ Az f z' "I L(6,z+Az) 

>i 1 T axr K ex P - sec e / z Kd " dS ' ~ T a M >i> 

i L 1 J 


z + Az 
i 


> 

r r z l 

T K exp -sec 0 j Kdu dz ‘ 
air r J 

z 

L i -J 


rZ + Az 
-sec 0J 1 Kdu 


sec 6i 


L{0 J z i +Az)f « K(e,co) 

= T b (0j v Az) - 


t cb 
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Therefore, 


T b <$, " T b ( 6 ’ z x + Az ) 


r 

air 


{z x ) - T b (0 jZi fA Z ) 


H 9, z + Az) 


or 


L{9, z + Az) 

146. FJ 


T b (0, Zi ) - T b {6, z^Az) 

T (z ) - T, (0, z + Az) 
air 1 ' b v 1 7 


(2-33) 


Hence, (2-32) and (2-33) form a pair of equations which may be used 

recursively to determine T (9, z+Az) and L,(G, z+Az) from the 

atm l i 

measurements when use is made of the known values of these quantities 
at the altitude z 

i 

For this procedure to be useful, an accurate calibration of the 
radiometers is required since (2-32) and (2-33) contain references to 
measured brightness temperatures which can be expected to be small 
under most circumstances. 

A second, more approximate, method of deteimimng atmospheric 
corrections may be useful when it is known that the atmospheric conditions 
of radiometric frequencies are such that the attenuation coefficient K is 

so small that sec 9 / Kdz «1. In this case, L>(0, z) 1 and 

J o 


T sky 


T 

atm 


(8,z) 


r m 

« sec fl K T dz + 

J Q air CB 

z 

« sec 0 / K T dz 
Jn air 


= sec 9 


1-7 

o o 


= hky- T CB- T b (e -“) 

This approximation can be used with data obtained m clear, low 
humidity weather 


(2-34) 
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Thus, only two measurements are required. The first, a measurement of 
I'sky near ground level, is analogous to the more piecise method 
discussed above. However, the measurements of the air temperature and 
brightness temperatures at a number of intermediate altitudes are 
eliminated and replaced by a single upward looking measurement of the 
brightness temperature, T-^(9, z) at the experiment altitude. As m the 
case of the more detailed procedure, the calibration of the radiometers 
must be sufficiently good so that a measurement of the brightness temperatures 
with the radiometers pointed skyward may be made more accurately 
than can be achieved with an a priori model. 
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Section 3 


ANTENNA PATTERN STUDY 

This section of the report documents equations and procedures 
developed to correct MFMR data for antenna pattern effects. 

3. 1 GENERAL CONSIDERATIONS 

An antenna detects energy propagating not only along its bore- 
sight axis hut also from other directions in space. The energy which 
is detected or, equivalently, the antenna temperature T^ may be 
expressed as 

T A ( e o>*o> = jG(e o ,cp o , e,<p)T r (G,<P) d a (3-1) 

where the angles ( 0^, cp ) specify the direction of the antenna bore sight 
axis with respect to some reference coordinate system, T r (9, cp) is the 
temperature of the radiation propagating along the ray specified by the 
angles (6, cp), dQ is a differential solid angle about (9, cp), and G is the 
antenna gain function. The gain function, which essentially describes 
the effectiveness of the antenna in detecting radiation propagaimg m 
directions other than the bore sight direction, may be normalized by 
requiring that 

Jg( 0 o , cp Q ; 9,cp)d0 = 1 (3-2) 


Since the cone of angles which contribute significantly to the 
antenna temperature will, in general, intersect the ground in an area 
which contains several different terrestrial features, it is clear that 
the antenna temperature will be a weighted average of many distinct 
ray temperatures. However, m the quantitative interpretation of 
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radiometric signals in terms of their relation to geophysical phenomena, 
it is the ray temperatures belonging to the distinct features which are 
desired instead of a weighted average. Thus, the problem of solving 
the integral equation (3-1) for the ray temperatures when the antenna 
temperature is known arises. 

There are some practical difficulties which severely limit 
any program designed to solve (3-1). First, there is the problem of 
measuring the antenna temperature at a sufficient number of points. 
Temperatures over a two-dimensional set are required ( 1 . e. , both along 
the aircraft flight path and perpendicular to the flight path). Since the 
antennas on the P3V are not scanned perpendicular to the flight path, the 
information which will be obtained will be confined to a narrow strip, 
which is essentially one dimensional. This is not adequate to solve (3-1) 
on the basis of microwave data alone. The situation can be improved 
somewhat by making several passes over the target area with the flight 
path displaced slightly on each pass. Another procedure for circumventing 
the lack of two-dimensional information is to make an a priori estimate 
of the temperature connected with features outside of the strip over which 
data is gathered. This may be done, for example, by an interpreter 
using optical data. 

A second problem that arises m connection with solving (3-1) 
and which is, perhaps, even more fundamental than the first is that 
of numerical instability. This problem will occur even if antenna 
temperatures are specified over a two-dimensional surface and will 
limit the degree to which the antenna beam may be "split" in any 
analytical scheme. The situation is as follows. In solving (3-1) 
numerically, the integral will be approximated by a finite sum of terms 
(obtained by the application of some numerical quadrature scheme) and 
the resultant system of linear equations with specified lefthand sides 
will be solved It will be found that if an attempt is made to resolve 
too many features by approximating the integral by a sum consisting of 
a large number of terms that a small change in the input data (i. e. , 
antenna temperature) will produce a large change in the computed ray 
temperatures. Thus, the exti action of meaningful ray temperatures 
becomes impossible. 
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The reasons for such behavior are common for all Fredholm 
integral equations of the first kind, of which (3-l) is an example m 
two dimensions. For simplicity in exposition, however, the problem 
will be discussed m one dimension m the following. Thus, we consider 
an equation of the form 

b 

T(x) = / G(x,y) T r (y)dy (3-3) 

a 


The limits of integration, a and b, are fixed. Suppose that the result 
T(x) is obtained with a given T_^(y). Now, let us examine the effect of 
modifying T^(y) by the addition of a term of the form A sin ny where 
A is a constant and n is a parameter. Then T{x) is changed by a 
term of the form 


A 


J* Gfe, y) 

a 


sin ny dy 


(3-4) 


However, in the limit where n -* oo, the integral in (3-4) approaches 
zero. This result follows rigorously from the Riemann Lebesgue 
lemma with certain mild restrictions on the class of functions from 
which G is taken (this class certainly includes functions which represent 
real antenna gain functions). Intuitively, the result is expected because 
the rapid oscillation of the term sin ny in the integrand will lead to 
cancellation between nearly equal positive and negative parts of the 
integral Thus, we see that even if A is a large constant ( 1 . e , there 
is a large change in T^(y) that the function T(x) is essentially 
unchanged when n is large. Convex sely, we see that a small change 
in T(x) (say, due to errors in measurement of the antenna temperature) 
maybe ascribed to large changes in the ray tempeiature. 


The argument presented above is not quite as artificial as 
it may seem at first sight in connection with ray temperature 
variations that will actually occur. It is only necessary to realize that 
the true temperature variation can be Fourier analyzed. Thus, we see 
that in the presence of small errors in measurement no information 
can be obtained about the high spatial frequency components of the ray 
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temperature. In fact, a numerical analysis using measured data 
can be expected to yield spurious, large high frequency components 
which will completely obscure the true nature of . Numerical 
experiments have indeed demonstrated that this phenomenon occurs 
and, unfortunately, is a serious problem even when only a small 
amount of "beam splitting" is attempted. Naturally, the problem 
manifests itself more severely as the error in the measured data is 
increased. 

In view of the above remarks, it is clear that stringent 
measures must be taken if meaningful antenna pattern corrections 
are to be made. This may be done by imposing certain smoothness 
requirements on the computed ray temperatures. Such a procedure 
can be used quite successfully if sufficient a priori information concerning 
the target characteristics is available. Unfortunately, this type of 
information cannot be expected to be available m general. Thus, another 
approach must be taken. Probably the simplest technique which can be 
used is to write (3-1) as 

T A (0 o> *<) = / GT r dfi + J GT r dn (3 

M S 

where M is the mam lobe of the antenna and S includes all of the 
sidelobes. A detailed description of the ray temperatures is not 
attempted. Instead, we attempt simply to obtain the average ray 
temperature in the mam lobe. The contributions of the sidelobes to 
the antenna temperature may be estimated by a procedure which is 
discussed in the next section. 

3. 2 DESCRIPTION OF CORRECTION TECHNIQUE 

The basic idea behind any sidelobe correction scheme Is to 
determine the value of the first integral m the expression 

T a - f T( 9, 0)G (9, 0)dO + J T(0, 0)G( 0, 0)dO (3 

M S 
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where is the observed antenna temperature, M represents the 

solid angle subtended by the mam lobe, S represents the solid 
angle subtended by the sidelobes, T(0, 0) is the brightness temperature 
at the surface m the direction (0, 0) and G{0, 0) is the gam of the 
antenna m the direction (0, 0). 

To solve the equation we must know G(0, 0) for each point 
m S and also T{^0) for each point m S where G{9, 0) is significant. 

In other words, we need to know the intensity of all significant sidelobes, and 
what these sidelobes were looking at 

The values of G(0, 0) for the MFMR antennas were determined 
by comparison with a gain standard 

Determination of T(0 f 0) for (0, 0) m S can be rather difficult, 
because the ground is scanned along a single line. Thus, on the basis of 
radiometric data, some correction maybe made for the sidelobes ob- 
serving the region immediately ahead and behind the mam lobe. However, 
it is unlikely that merely having data on what the sidelobes were seeing 
fore and aft of the mam beam will be sufficient to provide an acceptable 
correction. For most types of terrain, information will be needed on the 
terrain on all sides of the mam lobe It is expected that temperature 
values accurate enough for this purpose can be estimated using the optical 
photography obtained at the same time as the microwave measurements. 

Using the uncorrected microwave data and the photography, 
one would estimate the temperatures seen by the sidelobes within a 
specified angle of the ground track of the mam lobe. This would be done 
for whatever portion of the ground track that sidelobe -corrected data 
was required The estimates would then be entered onto punch cards 
or other computer -input media, and would be read into the computer 
along with the uncorrected microwave data. 
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The computer would then calculate the integral 


jT(6,0)G{e,0)dQ (3 

S 

for each point on the ground track, and subtract this value from the 
observed antenna temperature. The result would be the true main 
lobe temperature 


Jt( 0,0)G(4 0)dQ . 

M 

If these values differed excessively from the estimates used initially 
for a given type of terrain, a better set of estimates could be for- 
mulated and the procedure repeated. It is expected that after 
sufficient experience in preparing such estimates, the second 
iteration would not he required. 


(3 


3. 3 CORRECTION PROGRAMMING PROCEDURE 

The available information on the MFMR antennas indicates 
a mam beam efficiency of around 95% for the L-band (1.420 GHz) 
antenna, and around 90% for the higher frequency antennas. Either 
of these values is sufficiently below 100% to require a sidelobe 
correction procedure if the desired 1°K to 2°K absolute accuracy is 
to be achieved. 

It should he noted that all of the efficiency measurements 
were made by antenna range techniques, i. e. , by using a distant 
transmitter and receiving through the antenna under test. When 
this is done, the receiver is normally adjusted so as not to saturate 
on the peak main lobe signal. This means that sidelobe responses below 
about -40 dB will not be detected, and so will not contribute to the 
sidelobe portion of the mam lobe efficiency calculation. This would 
result m the main lobe efficiency being measured higher than it 
actually is This effect would be significant only if a substantial 
part of the far-out sidelobe pattern were above -30 dB. This is 
unlikely, but all flight data should be examined for any suggestion 
of such effects. 
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The available patterns seem to be adequate for determining the 
total power m the sidelobes. This is the case for the TRG antennas, 
because beam efficiency measurements made by TRG by use of ordinary 
measured patterns and symmetry assumptions agree fairly well with 
Aerojet measurements which were made by integrating the entire solid 
pattern with a pattern integrator. The 1 GHz antenna efficiency was 
determined by the pattern integrator method. 

3. 3. 1 COMPUTER PROGRAM OUTLINE 

a. The computer must be programmed to accept an approxi- 
mation of the microwave temperature of a strip of terrain centered on 
the ground track of the antenna beams and having a width equal to 
several times the aircraft altitude. Presumably, the ground track 
will either be straight or will curve sufficiently slow so that all 
points m the area of interest can be specified by a Cartesian cooidmate 
system with, for example, x representing distance perpendicular to 
the ground track and y representing distance along the ground track 
measured from an arbitrary reference - say, the beginning of the run 
Expressions for relating antenna beam position with ground position 
are derived m Section 3 4. 

b. The microwave temperature values for each point (x, y) 
would come from estimates prepared by someone familiar with the 
radiometers {and preferably familiar with the terrain) and having 
access to the uncorrected microwave data, as well as photographs, 
maps, etc. 

c. Since the antenna beams will, by definition, be tracking 
down the center of the specified terrain swath, the incidence angle 
will vary only with x, and then m a known manner. Therefore, the 
effects of varying incidence angle can be accounted for by specifying 
appropriately lower temperatures at larger values of x. 
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d. Significant portions of the sidelobe pattern will often 
point off the specified terrain region or even into the sky. For this 
reason it is necessary to specify "default options" giving an approxi- 
mation to the terrain temperature at shallow incidence angles xa 
various directions from the aircraft. The required number of 
different temperature versus incidence angle "curves" will 

depend on how variable the terrain is, but will probably not exceed 
4 or 5 for any given run. Also, the sky temperature will have to be 
specified as a function of incidence angle, and as a function of direction 
from the aircraft, when the sky conditions vary considerably with 
azimufh. The distance -terrain curves and the sky curves will vary 
from flight to flight, but can be selected (by the user) from a "catalog" 
of curves stored in the computer. 

e. In addition to specifying the radiometric nature of the 
terrain and sky it is also necessary to specify the time at which the 
antenna beams encounter a prescribed point on the Y-axis. This 
is done by specifying the time of encountering a small number of the 
points, and specifying the average ground speed between each pair of 
points. Also, the sidelobe pattern of each of the antennas must be 
specified and the incidence angle of the antenna beams fed into the 
computer. 

f. It will be most convenient to specify the sidelobe pattern 

in polar coordinates, centered on the main beam. Then the 6 coordinate 
will be "degrees off boresight" and the 0 cooidinate will he degrees 
circumferentially from some reference. When the antenna assembly 
is rotated to change polarization it is then only necessary to add 
90 ° to each 0 coordinate. The existing antenna patterns indicate 
that the cross-polarized response of all the antennas is negligible 
in the main lobe, and the aggregate response in the sidelobes can 
probably be accounted for with acceptable accuracy as a "bias" in 
the conventional (i. e., non-cross-polarized) sidelobe pattern. 

g. The parameter that would be used m specifying the 
sidelobe pattern is the fractional contribution to the total power 
received due to the given (0 , 0 ) "bin. " This depends on both the gain 
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at (0, 0) and the values of 0 and 0. The latter is necessary because 
the bin size, in steradians, varies with 0 and 0 while the gain is 
usually specified per steradian. However, it is only necessary to read 
in the usual gam pattern. The computer should be programmed to 
account for the size of each ( 0 , 0 ) bin. Also, since it is necessary 
to read the antenna pattern into the computer only once for each 
antenna, it is reasonable to use punch cards or something similar as 
the input medium for tins purpose Once the computer has the terrain 
information, the timing information, the incidence angle, the antenna 
pattern, and the microwave data available, the correction proceeds 
automatically as outlined m steps h through k and in the flow diagram. 
Figures 3-1 to 3-5. The flow diagram steps maybe followed exactly 
to generate the antenna pattern effects program. 

h. The computer considers each point in the antenna pattern, 

for example, starting with (just outside the main lobe) and 

incrementing 0 through 360°, then going to the next 8 "ring" and 
again incrementing 0 through 360°, and so on until 9 is reached, 

rrictx 

i. For each (0, 0) point where the gam is significant (say 
above -50 dB) the program determines where within the (x,y) terrain 
grid the (8, 0) direction intercepts the ground, or else winch one of 
the default -option terrain or sky temperature curves should be used. 

The program then could either use the temperature at that point (or 
nearest to that point) or, to obtain greater accuracy, the computer could 
determine the average temperature within the legion of the ground inter- 
cepted by the boundary of the (0, bin In either case, the temperature, 
multiplied by the fractional contribution of the given ( 0, 0) bin, yields 
the contribution to the observed temperature due to that bin. The 
computer adds the contributions from all the bins where the gam is 
significant, thus obtaining the total sidelobe contribution 

j. This total sidelobe contribution is then subtracted from the 
observed microwave temperature, and the result multiplied by the 
reciprocal of the fractional contribution of the main lobe to the total 
received power. The result is then the temperature seen by the main 
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INITIALIZE FREQUENCY 
FLAGS TO 1ST FREQUENCY 


INITIALIZE TIME TO 
MIN TIME 
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AT THIS TIME 


INITIALIZE S TO 0, (p TO 
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f CALCULATE 
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ST = sin (p CT = cos (f 

SX = SF*CT 
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SZ - (SA*SY-CF)/CA 
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.SZ < 0 9 
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SUBTRACT G($,4) (IN DB) 
FROM MAIN LOBE "GAIN 
ABOVE ISOTROPIC" (IN 
DB) FOR THIS ANTENNA, 
GIVING G a (0,0) 


3 

? 

CONVERT RESULT TO I 

NUMERICAL GAIN, J 

9 A (9r<t>) BY. 
g A (0,0) = 10 10 G A (0,0) 

% 

? 


CALCULATE VALUE OF d 
(ELEMENT OF SOLID 
ANGLE) FROM 

d42 = sin0d0 d 

(d0 AND dd ARE THE 
INCREMENTS FOR 0 AND 0) I 


3 ANTENNA PATTERN CORRECTION FLOW DIAGRAM 
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FIGURE 3-5. ANTENNA PATTERN CORRECTION FLOW DIAGRAM 
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lobe alone. 


k. This can be expressed mathematically as follows. 
The temperature observed at the antenna port is given by 



T M G M 


+ TjGj + T 2 G 2 + . . . 


(3-9) 


Where is the temperature seen by the main lobe, T^ is the 
temperature seen by bin 1 of the sidelobe pattern, T 2 is the temperature 
seen by bin 2, and so on. G^ is the fractional contribution to the 
total received power provided by the main lobe, Gj is the fractional 
contribution provided by bin 1, G 2 is the fractional contribution 
provided by bin 2, and so on. Solving for we get 

T M = 0- t A - < T 1 G 1 + T 2 °2 +•• ) 


The above procedure solves this equation exactly and will thus give 

the correct value of T^, provided the specified terrain and sky 

temperatures are correct. If, after running the progiam, it is 

found that the corrected main lobe temperatures for various terrain 

types fail to match the estimates used in the sidelobe correction 

procedure, then the newly-obtained information can be used to re- specify the 

terrain temperature and the program can be re-run. This procedure 

can be repeated as many times as necessary until the values converge. 

Presumably, after some familiarization with the radiometers and 

the terrain, only one or two runs will be required. 
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3.4 


DERIVATION OP ANTENNA BEAM POSITION EX- 
PRESSIONS 


The problem is to fmd the coordinates (x, y) on the ground 
representing the intercept of the direction (0, 0), where the (0, 0) 
coordinate system has its origin at (0, y Q jh) and its polar axis 
coincident with the bore sight axis (i. e. , in the YZ plane and elevated 
by an angle a from nadir toward the +Y direction). 

Let u^., Uy, and ti^ be unit vectors directed along the positive 
X, Y and Z axes, respectively. See Figure 3-6 

Let s^ be a unit vector directed away from the aircraft along the 
bore sight axis. 

Then 


-> 

s 

o 


sin cm - cos ou 
y z 


(3- 


Let s be a unit vector pointing from the aircraft toward the 
point (x, y, o) on the ground. 


T , — * -*■ . -> , -> 
Let s = su + su + su 
xx y y 2 z 


{Note s 2 + s 2 + s 2 = 1) 
v x y z ' 


Then, from the definition of dot products 


s«s 


= cos 6 = 


sm ors - cos os 

y z 


{3- 


Now, from the definition of cross products, the vector s > x s q is 

perpendicular to the plane defined by s and if , Also, s^ x (s”x if ) 

^ ^ o o o 

is perpendicular to both s and s x s . Therefore, s x (s x s ) is 

O ^ O O Q / 

m the plane defined by if and s q , and is also perpendicular to the boi*e- 
sight axis. 
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FIGURE 3-6. ANTENNA BEAM POSITION GEOMETRY 



From vector theorems* 


S 0 X ( SXS 0 ' “ ^ S 0 * S 0 ) S " ( S ’ S 0 ) S 0 = s ~ cos 0S Q 


Therefore, the unit vector zn the s^ x (s x s ) direction is* 


s - cos 0s 


s - cos 0s 


Now, 


s - cos 9s o = "y/(s - 


cos 0s )»(s - cos 0s ) 
O o 


=V;-; - 


2 cos 8s»s + cos 0s • s 
o o o 


= "Vl - COS^0 


Also, 


= sin 6 


COS 0 


s - cos 0s 


u • 

X 


s - cos 0 s 


u • s 

X 


sin 0 


sin 9 
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(3-13) 


(3-14) 


(3-15) 



Thus 


sin 9 cos 0 


(3-16) 


Now applying the fact that s^ + s^ + = 1 to eqation (3-11) 

we get: 


cos 8 = sin or s 

y 


From which* 


cos 8 - sin o' s 

y 


- cos or 



= -cos 


aV 


1 - s‘ 


x 


2 

s 

y 


Squaring both side s 

cos 9-2 sin or cos 0 s + sm 2 or s 2 = cos 2 or- cos 2 a s 2 - cos 2 cy s 

y y x ■ 

2 2 2 2 

cos 8 + s (sin or + cos 4 or) - 2 sin ffcos 8 s - cos 2 <y(l- s 2 ) = 0 
y y x 

But 

2 2 

sm or + cos or = 1. 


Thus, 

s 2 - 2 sin or cos 0 s + cos 2 8 - cos 2 cy (1-s 2 ) = 0 

y y x 

Solving the quadratic in s 2 * 

s y “ 1 2 sin or cos 0 ± 2 V sin 2 or cos 2 8 + cos 2 c<l- s 2 ) - cos 2 8 

X 


(3-17) 


(3-18) 
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Now by trigonometric identi.es 


* 2 2 Q 2 , , 2 \ 2 Q 2 Q 

sm of cos 0 - cos 0 = (1-cos o^cos 0 - cos 0 


2 2 

= -cos orcos 0 


Therefore 


2 2 2 2 

s - sin Qf cos 9 ± V-cos a cos 0 + cos a(l-s j 
y ^ x y 


2 2 2 2 
Taking cos or outside the radical and using the s = sin 0 cos 0 

we get 


Sy = sin O’ cos 8 ± cos cn ~V-cos^ 0 + 1-sm^ 9 cos^ 0 


- sm o' cos 0 ± cos o' sin 6 sin 0 

By continuity and by considering the situation near 0 = ir/2, we see 

that the + sign should be used. (This assumes that 0 runs from 0 to 2ir). 

Thus: 


s = sin a cos 8 + cos (x sin 9 sm 0 
Y 


and solving equation (3-11) for we get* 


sm of s - cos 9 

s = 1 

z cos a 


(3-19) 


(3-20) 


The vector s is now completely determined in terms of 0,0, and o'. 
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Now consider the line through the aircraft and the point (x, y, o), i. e. , 
the point where the (8,0) direction intercepts the ground. This 

) > t 

line is defined by the locus of points given by the vector + hu^ + 

where X is a parameter. Now the point on this line corresponding to 

(x, y, o) is the point for which the vector has zero z -component, i. e. , 

when X = -(h/s ). 

z 

•— > — > 

Now, the x and y components of y Q U^ + hu g + Xg are simply 
Xg^ and Xs^, respectively. Therefore, substituting -(h/s z ) for X we 
get the x and y coordinates of the desired point (x, y, o), namely 


hs 


x - 


x 


(3-21) 


hs 


+ y. 


(3-22) 


The criterion for the (8,0) direction intercepting the ground is obviously 
s_ <0. If this is not met, or if the x or y coordinates are beyond 
the area for which temperatures have been predicted, then we can use 
either zenith angle and aximuth to find the sky temperature m the (0, 0) 
direction or depression angle and azimuth as a guide to finding the 
brightness temperature of the ground 


The zenith angle is simply the angle between s and the + z -axis, 
which is cos * s . The depression angle is ir/2 less than the 
zenith angle, namely, (cos' s^) - tt/ 2 and the azimuth is tan {-s /s x ) 
measured counterclockwise from the Y-axis, (i.e., counterclockwise 
from the direction of the nose of the aircraft). 


It should be noted that the FORTRAN function 
equivalent) should be used to calculate tan -1 (-Sy/s x ) 
the angle is in the correct quadrant. 


ATAN2 (or 
to insure that 
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Sectipn 4 


EXPERIMENTAL VERIFICATION OF CORRECTION ROUTINES 
4. 1 GENERAL 

This section documents experimental tests developed to verify the 
adequacy of atmospheric and antenna pattern correction routines outlined m 
Sections 2 and 3 of this report. These experiments included bench and in- 
flight calibration of the MFMR system, along with airborne measurements 
to evaluate the correction equations. 

MFMR missions were flown on 13 and 16 February 1970 m an 
attempt to verify the correction routines. These flights were over the Gulf 
of Mexico and the Houston area Cold and hot load calibration data acquired 
during late January 197° and late February 1970 were utilized m conjunction 
with mission data. Experimental data taken to substantiate the atmospheric 
correction routines show that gross inconsistencies exist m the MFMR data. 
Unrealistic brightness temperature values were obtained using the antenna 
and radome loss values furnished for this study Significant improvement was 
achieved by using an empirically determined loss factor. However, quanti- 
tative results could not be obtained due to the sensitivity of measurement 
data to the empirically determined loss factor This uncertainty, discussed 
in detail m Section 4. 2, must be resolved before the reported atmospheric 
and antenna pattern routines can be substantiated. 

4. 2 EXPERIMENTAL VERIFICATION OF ATMOSPHERIC 

CORRECTION ROUTINE 

4. 2. 1 MISSION PLAN 

An experimental mission was configured to test the validity of 
atmospheric correction routines developed m Section 2 The mission was 
approved by cogmzant NASA/MSC personnel and is outlined below. 
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ATMOSPHERIC MISSION PLAN 


A. Flight Conditions 


1) Flight should be over ocean, prefereably having constant 

sea state characteristics. A cloudless day with calm seas is 
preferred. 

2) Flights should be level and repeated passes made over the 
same flight line, in the same direction for stated angle 
scans. 

3) Flight line one (1) consists of repeated passes as shown on 
Table 1. Total passes = 40. 

4) Radiometers pointing downward (0°) except for part of low 
level run and part of high level run. See Table 4-1. 

5) Record vertical polarization only. 

6) Record output during 500 ft elevation and ~ 20, 000 ft 
elevation pass when swinging from vertically downward (0 ) 
to ~ vertically upward (150°) m 10° steps. Record for 

30 seconds at each 10 u increment. 

7) Record baseline and cal "1" before and after each flight 
run. (10 secs each.) Note The cal "after" a given run 
becomes the cal "before" the next run, etc. 

8) Bore sight camera imagery is required during each data run. 

9) Repeat all steps in a descending flight 

10) Estimated total flight time, both ascending and descending - 
80 to 100 minutes. 

11) Low level pass, 500 ft, high level pass, 20, 000 ft; steps 
(ascending), approximately 1100 ft, steps (descending), 
approximately 3, 000 ft* 

12) Elevation control should be by radar and not barometric 
pressure. 

13) Integration times should be 1.0 sec for 10. 625 GHz 
radiometer (X-Band) and 0. 1 sec for the other radiometers. 
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ATMOSPHERIC MISSION PLAN (FLIGHT LINE 1) 


TabJ« 4-1 

ATMOSPHERIC MISSION PHAN (Flight Line 1) 
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Run 

No 

Flight 

Line 

Pass 

Absolute 

Flight 

Altitude*** 

Record 

Time 

1 

1- 1 

500 

30 i 

secs 

2 


500 

10 mins 

3 


500 

30 i 

9ecs 

4 

1- 2 

1, 600 


fc 

5 

1- 3 

2, 700 



6 

1- 4 

3 800 



7 

1- 5 

4,900 



8 

1- 6 

6, 000 



9 

1- 7 

7 100 



10 

1- 8 

8, 200 



11 

1- 9 

9, 300 



12 

1-10 

10 400 



13 

1-11 

11 500 



14 

1-12 

12, 600 



15 

1-13 

13 700 



16 

1-14 

14 800 



17 

1-15 

15,900 



18 

1-16 

1 700 x 



19 

1-17 

18, 100 



20 

1-18 

19, 200 

t 

f 

21 

1-19 

20, 300 

l 

30 i 

IOCS 

22 


20 300 

10 mins 

23 


20, 300 

30 secs 

24 

1-20 

19, 200 


i 

25 

1-21 

18, 100 



26 

1-22 

17 000 



27 

1-23 

IS 900 



28 

1-24 

14, 800 



29 

1-25 

13,700 



3D 

1-26 

12 600 



31 

1-27 

11,500 

30 secs 

32 

1-28 

10 400 



33 

1-29 

9 300 



34 

1-30 

8 200 



35 

1-31 

7 100 



36 

1-32 

6, 000 



37 

1-33 

4 900 



33 

1-34 

3, 800 



39 

1-35 

2 700 



40 

1-36 

1, 600 


7 

41 

1-37 

500 

30 secs 

42 


500 

10 mins 

43 


500 

30 secs 

Note * 

Sensors include multfrequency radic 


Sensors* 

Yes 

A 


Laser 

Profilo- 

metar 

Yea 

A 


Meteoro- 

logical 

Data 

Yes 

A 


Polarl- 

zation 


? 

? 

0 

A 


Yes 


V 

yes 


V 


V 

yea 


Vert 

A 


Yes 


Vert 


V 

yea 


** 0 view angle is vertically downward, 150° view angle is towards sky 

*** Altitudes may approximate ±200 given altitudes however exact altitude of pass 
must be recorded 


View Angle 

0° ** 

0° to 150 in 
10° Steps 

150° 

0 ° 

A 


* 

0 ° 

0° to 150° in 
10° steps 

150° 

0 ° 

A 


V 

o° 

0° to 150° 
in 10° steps 
ISO® 
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B. Recorded Flight Data 


1) Radiometer outputs and radiometer calibration data on 
PCM system and on 5-pen analog recorder (use 1. 0 sec 
integrators) (record 4 frequencies, plus time). 

2) Radiometer Frequencies 


L-Band 

1.42 

GHz 

X-Band 

10. 625 

GHz 

K-Band 

22. 2 

GHz 

K -Band 

31. 4 

GHz 


3) Complete ADAS 

4) Bore sight camera imagery during all recording periods. 
(Must be time or frame correlated with radiometer data.) 

5) Air temperature at each elevation 

6) Air pressure at each elevation 

7) Humidity at each elevation 

8) Occurrence of water droplets (if any) at each elevation. 
Quantitative cloud cover. Note height of cloud base and 
cloud tops in flight line (if any). 

9) Record laser altimeter data at applicable lower elevations. 
C. Ground Truth Data 

1) Nominal wind speed (steady or gusting) 

2) Nominal wind direction 

3) Nominal wave height 

4) Nominal water temperature 

5) Description of weather m general area. (Meteorological 
maps ) Presence of fronts, cloud cover descriptions, 
ceiling (if any). 

6) Cold load and hot load calibrations should be run within a 
day or two "before" and "after" the flight. 
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4.2 2 


FLIGHT TEST RESULTS 


A mission was flown on February 13, 1970 m an attempt to verify 
the atmospheric correction equations described in Section 2.2 1. Data were 
gathered over the Gulf of Mexico with the radiometers pointing both toward 
the sea and skyward. Of particular interest from the point of view of testing 
the basic equations without the extraneous influence of emissions from other 
sources is the case where the radiometers were pointing skyward. These 
data also provide a good test of the radiometer calibration (including radome 
effects) since low temperature targets generally provide sensitive indications 
of possible calibration errors. The following discussion will be restricted to 
a consideration of only this portion of the flight test data since the gross 
inconsistencies which were encountered seem to preclude a precision esti- 
mate of temperatures from other targets until corrective steps are taken. 

The sky temperature measurements can be conveniently divided 
into two groups for the purpose of discussion. The first part of the data was 
gathered as the aircraft climbed from an altitude of 550 to 21, 000 feet, while 
the second part of the data was obtained during a descent from an altitude of 
19, 900 to 485 feet. In the following, these portions of the flight will be 
referred to simply as the ascent and descent. In both cases the air tempera- 
ture and dew point (from which the absolute humidity can be determined) 
were recorded at a number of altitudes. No clouds weie present. 

The measured air temperatures for the ascent and descent are 
shown m Figure 4-1. Since the calculation of sky brightness temperatures 
requires atmospheric data at altitudes higher than were attained by the air- 
craft, an extrapolation is required, A glance at Figure 4-1 shows that the 
measured air temperatures below 6 km do not differ gi eatly from the ARDC 
model tropical atmosphere. Hence, m the calculation of sky brightness tem- 
peratures, measured air temperatures were used from sea level to the peak 
flight altitude (slightly greater than 6 km), between the peak flight altitude 
and 8 km, a smooth transition into a tropical air temperature profile was 
introduced, above 8 km, the model tropical atmosphere was used. 
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The absolute humidity, p H Q (z) , during the ascent and descent is 
shown m Figure 4-2. A few of the measured dew points lead to relative 
humidities which are greater than 100 percent {assuming that the recorded 
air temperature is correct). There is clearly some error here and the 
points which are plotted (indicated by question marks m the figure) are the 
absolute humidities assuming 100 percent relative humidity. Both the ascent 
and descent measurements indicate that the water vapor density decreases 
approximately man exponential manner with altitude. Thus, since informa- 
tion at alitutdes higher than were achieved in the flight is needed for a theo- 
retical calculation, the assumption was made that the water vapor profile is 
adequately described by 


p h 2 o ( z ) ~ p h 2 o (0) 6 


For both the ascent and descent, p^. q( 0) was chosen to be 22.0 gm/mC 
while a reasonable fit to the data is given by \ = 0.523 km * fox the ascent 
and A. =0.618 km ^ for the descent portion of the flight. 


On the basis of the model discussed above, sky bughtness temper- 
atures were calculated for frequencies of 1,42, 10.625, 22.235, and 31.4 
GHz at a number of viewing angles and aircraft altitudes. The results are 
shown m Figures 4-3 to 4-12. The uncertainty m the calculated values due 
to the required extrapolation of air temperature piofile is small. There is 
more uncertainty due to the approximate water vapor density profile which 
was used. The uncertainty is greatest at 22.235 GHz (near the peak water 
vapor absorption) and negligible at 1.42 GHz Except near grazing (i.e., 
angles less than 110°), the uncertainty at 22.235 GHz should be less than 
10°K near sea level and less than 2°K near 20, 000 feet. 


These results may be compared with measured brightness tempei - 
atures. Comparisons, using several different procedures m processing the 
data, are shown m Figures 4-13 to 4-15 for frequencies of 1.42, 22.235, 
and 31.4 GHz. 
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FIGURE 4-2 WATER VAPOR DENSITY PROFILE - ATMOSPHERIC EFFECTS 
TEST FLIGHT, 2/13/70 
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FIGURE 4-4 T b VS. VIEW ANGLE AND ELEVATION 
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T b vs VIEW ANGLE AND ELEVATION 
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In all of the discussion which follows we make the assumption that 
no significant error exists m the measured radiometer data. The validity of 
this assumption can only be determined after extensive measurements have 
determined the errors m each parameter used to obt am brightness temperature 

The first method of deriving results from the raw data consisted of 
applying the calibration data for the radome and antenna which were obtained 
at Table Mountain^. See Section 4. 2. 3. At all frequencies, totally unrea- 
listic results were obtained. The most obvious errors are the large negative 
temperatures at 22.235, 22 355 and 31.4 GHz. Negative temperatures are 
physically meaningless. 

Another difficulty that occurs with the Table Mountain calibration 
data is the oscillation m the sky temperature with angle (recall that we have 
assumed the measured data to be correct). In the absence of a region of 
strong absorption at a particular angle, the sky temperature should be a 
raonotomcally increasing function of angle m scanning from zenith to grazing (l 
from 180° to 90°). This behavior was not obtained m the measurements and 
cannot be ascribed to a possible local atmospheric anomaly because the air- 
craft was heading m different directions at different times, thus allowing the 
radiometers to observe several different portions of the atmosphere. Finally, 
the use of the Table Mountain calibration leads to an unexplamably high bright- 
ness temperature (above 70°K) at 1.42 GHz. 

In view of the failure to obtain meaningful results with the pre- 
scribed calibration data, a second, approximate method was tried m order 
to test the sky temperature model. This consisted of empirically determin- 
ing an overall loss factor (see Section 4.2. 3) at each frequency and at an 
angle of 150° by forcing the sky bughtness temperatures at 21, 000 feet to 
agree with the theoretical results. This loss factor was then used to reduce 
the data at other angles at 21, 000 feet as well as at 550 feet (for the ascent) 
and 485 feet (for the descent). 


(1) "Table Mountain Calibration on P3A Multifrequency Microwave Radiom- 
eters and Radomes, " Jet Propulsion .Laboratory, Nov. 25, 1969. 
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With the second method of calibration, there is a distinct improve- 
ment m agreement with theory at all frequencies. The measurements near 
sea level are of the correct order of magnitude. Quantitative agreement, of 
course, cannot be expected because of the sensitivity of the results to the 
empirically determined loss factor. This factor could not be fixed precisely. 
The principal qualitative failure of the second calibration method is the pecu- 
liar oscillation m the brightness temperature which was previously noted m 
the use of the first calibration method. It must be concluded that allowance 
must be made for a loss factor which varies with angle. 

Thus, a third approach was tried. This consisted of determining a 
loss factor as in the second method with the exception that a separate loss 
factor was assigned to each angle With this method, all of the data at 
El, 000 feet were used to establish the calibration so that the only check with 
theory is the data at lower altitudes ( 1 . e. , the 550 foot ascent data and the 
485 foot descent data). The results are shown m Figures 4-16 to 4-18. A 
substantial improvement is made over the second method m that the magni- 
tude of the undesirable oscillations is greatly decreased. Again, because of 
the sensitivity to the loss factor and approximations which were necessary, 
a precise quantitative check cannot be expected. 

4.2.3 DATA REDUCTION PROCEDURES 

In the computation of brightness temperatures from data recorded 
by the data acquisition system used m the flight test measurements, a num- 
ber of operations take place. This section is concerned with presenting, m 
some detail, the operations which aie needed to arrive at radiometer tem- 
peratures after the measured data has been supplied to Aerojet -General by 
the Computation and Analysis Division of the Manned Spacecraft Center m 
Houston, Texas. 

The data acquisition system used m the aircraft samples and 
records each parameter of interest approximately 60 times per second. The 
values of each parameter can then be averaged using a computer program 
available at the Spacecraft Center with the time of measurement serving as 
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Legend, Ascent 
O 500 ft. Theoretical 
□ 21. 1 K ft. Theoretical 
V 500 ft, Lata 
A 2 1 1 K ft. Data 
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FIGURE 4-16 

T b VS VIEWING ANGLE, VARIABLE LOSSES 
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a means of identifying the data. Before the computer program can perform 
an average,however, two quantities must be supplied a time interval which 
specifies the number of values to be averaged and the time rate of printing 
the averaged values. The average value is printed out m the mid -point of the 
time interval. For example, when a five-second time interval and one-second 
time rate of printing is employed, 300 successive values of a parameter are 
used m the average with the average value printed out every second at the 
time corresponding to the mid -point of the five-second interval. For the data 
used m this report, a five-second time interval and a five-second time rate 
of printing was employed. 

To readily determine the location xn time of possible large noise 
fluctuations in the recoided data, an analog output on chart paper of the 
recorded data was examined. Locations which indicated large noise fluctua- 
tions were subsequently not considered as representing meaningful data. 


Using averaged values of measui ed parameters supplied by the 
Manned Spacecraft Center, uncorrected brightness temperature T^ can be 
calculated as a function of time t 


T b<*> = 


T nt( V R< t > - W’) 
V c<*> - V BL<« 


+ T ref<« 


where corrections for antenna and radome losses occur later 


and 


In the above equation 
T^ 1S the temperature of the noise tube, 

T^ gf is the reference tempeiature of the baseline (°K), 
V c is the calibrate voltage, 

V-gk 1S the baseline voltage, 

is the radiometer voltage. 
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Measurements of the temperature of the noise tube were made 
before and after the flight test measurements by O'Hare'*' for all frequencies 
with the exception of pre-flight measurements at 22.235 and 22.355 GHz. 
Table 4-2 gives the values m °K As can be seen from the table, a rela- 
tively large uncertainty m exists for 1.420 and 10.625 GHz, 

Table 4-2 


Frequency 

(GHz) 

Pre -Flight 
Measurements 

Post -Flight 
Measurements 

Average 

1.420 

123. 1 

119.4 

121.2 

10.625 

94,4 

89.6 

92.0 

22.235 


107.3 

107.3 

22.355 


124.2 

124.2 

31.400 

85.8 

85.5 

85.6 ; 

i 


The reference temperatures of the baseline during the flight meas- 
urements are shown m Figures 4-19 through 4-22 for each frequency. The 
error m the temper atui es is usually less than 0,2°C for the theimistors 
employed (YSI precision thermistors weie used). Variations greater than 
0. 1°C of the reference temperature dunng the flight measurements were 
included m the calculations of T^ . 

The baseline and calibrate voltages were monitored before and 
after each radiometric measurement The maximum and minimum values at 
a fixed point in time of the five-second averaged data was observed to be 
approximately 30 millivolts for the 1.420 GHz radiometer. Variations of the 
same order of magnitude also exist at the other frequencies. 


*> Interdepartmental communication from W, M, O'Hare to N. M. Hatcher 
concerning MFMR calibration data on 3/30/70 at Manned Spacecraft 
Center, Houston, Texas. 
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FIGURE 4-20 REFERENCE THERMISTOR (T^) VOLTAGE VS TIME 
10 625 GHz (X-Band), 2/13/70 


rnermibtor Voltage (Volts; 



FIGURE 4-21 REFERENCE THERMISTOR (T^ VOLTAGE VS TIME 
22 235 and 22 355 GHz (K p K z Band), 2/13/70 




1 


FIGURE 4-22 REFERENCE THERMISTOR (T x ) VOLTAGE VS TIME 
31.4 GHz (K Band), 2/13/70 

a 


Since the variation of the average value of and during the 

flight were roughly the same as the variations occurring at a fixed point m 
time, the variations of the average during the flight were considered to be due 
noise fluctuations rathe r than as real fluctuations of and . Straight 

lines were constructed to pass through the average values and used m the 
calculation of T^ . See Figures 4-23 through 4-27. Each point m the figures 
is the average of more than 900 sampled values. It should be noted that, at 
times, values of calibrated voltages were not recorded by the PCM data sys- 
tem m the flights. In these cases linear extrapolations were used to approxi- 
mate the calibrate voltages. 

From an observation of the variations of V and Vt-, t m the 

c BL 

figures it is apparent that detailed long-term measurements of these quanti- 
ties would certainly be informative. The measurements would determine 
what portions of long-term variation of measured values of V c and 

are due to real variations of V and V„ T . 

c BE 

Table 4-3 gives values of T^Ct) corresponding to average values 
of V^(t) together with an rms error due to fluctuations m . The num- 
ber of points used m determining the average and rms is also given. As seen, 
the rms error is relatively large for tie 22 235 and 22.355 GHz radiometers. 
Note also that two inexplicable large rms values occui at 31.4 GHz. As 
observed m the table, the rms error at 10.625 GHz is often 0.0. A prelim- 
inary observation indicates the possibility that voltages less than zero out of 
the radiometer are somehow recorded as zero volts. This possibility should 
be investigated. 

Corrections of brightness temperatures for effects of losses of 
radome and antenna are usually made using the equation 


T 

correct 


(t) 



T b (t > - W> - V 1 - l a> t r<« 


where and are the power loss fractions of the radome and antenna 

respectively and T-^ and T^ ai e the thermal temperatures of the radome 
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Table 4-3 


UNCORRECTED BRIGHTNESS TEMPERATURES 
AT THREE ALTITUDES* 


Elevation 
& Angle 

1.4200 GHz 

10.625 GHz 

22 235 GHz 

22.355 GHz 

31.400 GHz 

No. of 
Points 

Ascent 
550' & 100° 

172.3 ±0.9 

59.9 ±0.4 

130.2 ±1.7 

153.2 ±1.0 

151.3 ±0.4 

8 

550’ & 110° 

164.7 ±0.9 

53. 3 ±0. 1 

84.6 ±1.6 

110.5 ±1.7 

141.7 ±0.2 

7 

550' &120° 

165.2 ±0.5 

53 7 ±0. 1 

77.8 ±1.4 

104 8 ±1.4 

122.8 ±30.5 

7 

550' & 130° 

164.1 ±0.5 

53.8 ±0. 1 

61.8 ±1.6 

91.5 ±1.5 

134.8 ±0.7 

5 

550’ & 140° 

163. 1 ±0.9 

54. 1 ±0.1 

68.8 ±1 7 

96.4 ±1.3 

116.0 ±0.7 

6 

550* & 150° 

163.0 ±0.5 

54.6 ±0.3 

76.2 ±1.9 

101.7 ±1.0 

113.5 ±0.4 

5 

, 100’ & 100° 

157.5 ±1. 1 

65.5 ±0.4 

15.6 ±1.6 

51.5 =0.7 

96.8 ±0.4 

9 

, 100’ L 110° 

153.0 ±0 4 

61.6 ±0.0 

8.9 ±1. 1 

47. 3 ±1.0 

103.2 ±0.3 

4 

, 100' & 120° 

152.6 ±0.8 

61.8 ±0 0 

22.7 ±1.4 

57.9 ±2. 1 

99.5 ±1.0 

4 

, 100' &130° 

151.4 ±0.6 

61.9 ±0 0 

20.6 ±1.8 

53.7 ±1.0 

99.8 ±0.4 

4 

, 100' & 140° 

151.2 ±0.5 

62.0 ±0.0 

28. 1 ±1.2 

62.4 ±0.8 

71.9 ±16.0 

3 

, 100' & 150° 

149.7 ±1.0 

62.0 ±0.0 

39.5 ±1.0 

71.6 ±1 1 

84.8 ±0.5 

3 

Descent 
485' & 100° 

171. 1 ±0.3 

69. 6 ±2. 5 

154.4 ±7.7 

177.5 ±1. 1 

159.5 ±0.5 

6 

485' & 110° 

162.6 ±0.7 

69.0 ±1.1 

108.6 ±1.5 

133.2 ±1.2 

147.2 ±0.5 

5 

485’ & 120° 

160.6 ±0.4 

68.8 ±0.0 

98.2 ±1 0 

123.3 ±0.7 

139.7 ±1. 1 

5 

485’ & 130° 

159.9 ±0.5 

68 8 ±0.0 

85.9 ±1.8 

110.5 ±1.2 

137. 1 ±0.3 

4 

485' & 140° 

161.2 ±0 6 

69.2 ±0.0 

94 7 ±1 4 

118.3 ±1.9 

119. 1 ±0.2 

4 

485' & 150° 

162.5 ±0.7 

69. 3 ±0.0 

101.5 ±1.3 

124.6 ±1.2 

118.6 ±0.7 

5 


* The rms error is calculated from fluctuations m . The number of sample 
values used m arriving at an average V R is 300 times the number of points 
given m the table. 
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and antenna respectively. A number of assumptions are made m using this 
relation 

a. Reflections from radome or antenna are negligible. 

b. Thermal gradients m the radome and antenna are negligible. 

c. Variations of losses mthe radome and along the antenna 
are negligible. 

An investigation of the validity of these assumptions would allow 
an estimate of the error introduced into brightness temperatures using the 
above equation. 

The values of and L^ were taken from the Table Mountain 
calibration and are listed below*. 


Table 4-4 

LOSSES OF RADOME AND ANTENNAS 


Frequency 

L^ (Radome) 

L^ (Antenna) 

1.420 

0.077 

0.299 

10.625 

0. 107 

0. 160 

22.235 

0.286 

0. 190 

22.355 

0.268 

0.272 

31.400 

0.269 

0. 168 


The thermal temperature of the radome was monitored by six 
thermistors located on top outside, left outside, left middle, left inside, 
bottom outside and bottom inside of the radome. Table 4-5 gives the average 
values measured at two elevations. An average of the readings at the corre- 
sponding altitude was used mthe correction of measured brightness temperatur 

4 Table Mountain Calibration on P3A Multifrequency Microwave Radiom- 
eters and Radomes, Nov. 25, 1969, JPL. 
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Table 4-5 


THERMAL TEMPERATURES AT DIFFERENT 
LOCATIONS AROUND THE RADOME IN °C 


Elevation 

Top 

Outside 

Left 

Outside 

Left 

Middle 

Left 

Inside 

Bottom 

Outside 

Bottom 

Inside 

Sea Level 

23.4 

24.4 

23.7 

23. 1 

21.0 

21.4 

21, 100 feet 

-0.8 

0 
• 

ID 

1 

-2.8 

-2.4 

-5.8 

-4.0 


The thermal temperature of the antennas and waveguides were 
monitored with two thermistors. Although RF loss occurs m the waveguide 
and coaxial leads, most of the loss occurs m the antennas. It is for this 
reason that the thermal temperature of the antenna is used with , The 
thermal temperature of the waveguide is usually only a few degrees warmer 
than the antenna so that negligible error is introduced m using the thermal 
temperature of the antenna 

Table 4-6 below gives the corrected brightness temperature cor- 
responding to the uncorrected values m Table 4-3. As mentioned m the dis- 
cussion of results, gross discrepancies occur between the theoretically pre- 
dicted values and those m Table 4-6. 

If one assumes that the measured data has no significant error, 
then one can calculate an equivalent loss of the entire radome-radiometer 
system. That is, if the quantity L^ 1S selected from 



where T^ is the measured brightness temperature of the sky, T g is the 
theoretically predicted value and T' is some average thermal temperature 
of the radome and antenna, then all other measured data can be corrected 
for losses using 
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Table 4-6 


CORRECTED BRIGHTNESS TEMPERATURE 
USING THE LOSS VALUES DETERMINED ON TABLE MOUNTAIN 
(See Table 4-4 for values of L^ and V 


Elevation 
& Angle 

1 4200 GHz 

10.625 GHz 

22 235 GHz 

22.355 GHz 

31 400 GHz 

Ascent 
550' & 100° 

104.5 

-17.9 

10.2 

21.8 

57.8 

550' & 110° 

92.8 

-26.7 

-64.6 

-36. 6 

42. 1 

550' & 120° 

93.6 

-26 2 

-76.0 

-44. 3 

31.7 

550- & 130° 

92.0 

-26.0 

-102.0 

- 62.6 

30.9 

550' & 140° 

90.4 

-25 6 

-90.6 

-55.8 

0.0 

550' & 150° 

90.2 

-25.0 

-78.6 

-48.6 

-3.9 

21, 100' & 100° 

94.5 

16.8 

-173.0 

-104.6 

-17.4 

21, 100' & 110° 

87.6 

17 0 

-185,0 

-110.5 

-6.9 

21, 100’ & 120° 

86.7 

17 2 

-161.0 

-96.0 

-13.0 

21, 100' & 130° 

85. 0 

17.4 

-164.4 

-101.8 

-12.5 

21, 100’ & 140° 

84.7 

17.5 

-151.5 

-90.0 

1 

o 

• 

Descent 






485' & 100° 

104.0 

-4.4 

52.7 

55.0 

73.4 

485' & 110° 

91.0 

-5.2 

-26.4 

-5.3 

54. 1 

485' & 120° 

88.0 

-5.5 

-44.5 

-18.8 

41.8 

485' & 130° 

86.7 

-5.5 

-65.5 

-36. 3 

37.6 

485' & 140° 

88.8 

-4 9 

-50.5 

-25.6 

8.2 

485' & 150° 

90.8 

-4.8 

-38.8 

-17.1 

7.2 
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correct 


1 - L 


T, 


eff 


L ,, T' 
eff 


The theoretically predicted sky temperature is computed from 
measurements of water vapor content, water droplet content, relative 
humidity, pressure and thermometric temperature profiles of the atmos- 
phere. See Section 4. 2. 2. 

In the figures comparing measured to corrected brightness tem- 
peratures with theoretically calculated temperature, the notation "sky 
adjust losses" refers to temperatures calculated using L^ • 

Values of L^ calculated from sky measurements taken at 
21, 100 feet are given m Table 4-7 where T ' was determined using an aver- 
age of the radome and antenna temperatures. The values were determined 
using the sky measurement at 150° from nadir. 


Table 4-7 

VALUES OF L „ DETERMINED AT TWO ELEVATIONS * 
eff 


Frequency, GHz 

L „ 
eff 

550 Feet Elev 

485 F eet Elev 

1.420 

0.542 

0.542 

10.625 

0.215 

0 216 

22.235 

0. 104 

0. 118 

22.355 

0.226 

0.238 

31.400 

0.292 

0.295 


* The 550 foot measurements corresponds to ascent data while the 485 foot 
to the descent data. Also, note that L g ^ could be determined at each 
angle of observation. This was done and the results are discussed in 
Section 4.2.2. 
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4.3 


EXPERIMENTAL VERIFICATION OF ANTENNA PATTERN 
EFFECTS CORRECTIONS 


An antenna pattern effects experiment was designed to provide 
information suitable for verifying the techniques of the antenna pattern cor- 
rection routine outlined in Section 3 and for demonstrating the magnitude of 
such corrections on brightness temperature data. The test design is described 
m Section 4. 3. 1 and flight testing was conducted on 16 February 1970. The 
antenna pattern flight test data were not analyzed because of the sizable 
uncertainties m MFMR data discussed m the preceding section These 
uncertainties preclude a meaningful analysis of the flight test data. 


4. 3. 1 MISSION PLAN 

An experimental mission was configured to test the validity of 
atmospheric correction routines developed m Section 3. The mission plan 
was approved by cognizant NASA/MSC personnel and is outlined below. 


ANTENNA PATTERN MISSION PLAN 
A. Flight Conditions 

1) Flight should be over areas having sharp demarcation of 
media. Viz. water/land, land/concrete. Pairs of passes 
{using horizontal and vertical polarization) should be made 
over same flight line at six different altitudes. These 
flights should be accomplished using an antenna viewing 
angle of 0°, 30° and 150°, as indicated m Table 4-8. 

2) Flight should be on a clear day. Time of day is unimportant. 

3) Flights should be level and repeated passes made over the 
same flight line, m the same direction. 

4) Flight passes should be perpendiculai to boundaries. 

Record approxomately 30 seconds before and after each 
boundary interface. 

5) Synchronize boresight camera imagery with each data run. 

6) Specific flight lines will be m Houston area (within radius of 
250 miles) and will be precisely defined prior to 12/31/69. 
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TabLe 4-8 


ANTENNA PATTERN MISSION PLAN 
(Flight Lines 2 and 3) 


Run 

No. 

Flight 

Line 

Pass 

Absolute 

Flight 

Altitude 

Record 

Time 

Sensors 

1 

Laser 

Profi- 

lometer 

Meteoro- 

logical 

Data 

Polari- 

zation 

Ant enna 
View 
Angle 

1 

2- 1 

250 

60 secs 

Yes 

Yes 

Yes 

Hor. 

0° 

2 

2- 2 

250 

A 

k 



A 


A 


Vert. 

0° 

3 

2- 3 

250 









Hor. 

30 o 

4 

2- 4 

250 









V ert. 

30° 

5 

2- 5 

250 









Hor. 

150° 

6 

2- 6 

250 









Vert. 

150° 

7 

2- 7 

500 









Vert. 

150° 

8 

2- 8 

500 









Hor. 

150° 

9 

2- 9 

500 









Vert. 

30° 

10 

2-10 

500 









Hor. 

30 o 

11 

2-11 

500 









Vert. 

0° 

12 

2-12 

500 









Hor. 

0° 

13 

2-13 

1000 









Hor. 

0° 

14 

2-14 

1000 









Vert, 

0° 

15 

2-15 

1000 









Hor. 

30° 

16 

2-16 

1000 









V ert. 

30° 

17 

2-17 

1000 









Hor. 

150° 

18 

2-18 

1000 









Vert. 

150° 

1° 

2-19 

2000 









Vert. 

150° 

20 

2-20 

2000 









Hor. 

150° 

21 

2-21 

2000 









V ert. 

30° 

22 

2-22 

2000 









Hor. 

30° 

23 

2-23 

2000 









Vert. 

0° 

24 

2-24 

2000 









Hor. 

0° 

25 

2-25 

4000 









Hor. 

°o 


2-26 

4000 









Vert. 

0° 

27 

2-27 

4000 









Hor. 

30° 

28 

2-28 

4000 





4 

7 



V ert. 

30° 

29 

2-29 

4000 





Yes 



Hor. 

150° 

30 

2-30 

4000 








V ert. 

150° 

31 

2-31 

8000 








V ert. 

150° 

32 

2-32 

8000 








Hor. 

150° 

33 

2-33 

8000 








Vert. 

30° 

34 

2-34 

8000 








Hor. 

30° 

35 

2-35 

8000 


7 

* 

f 


4 

f 

V ert. 

0° 

36 

2-36 

8000 


Ye 

5S 


Yes 

Hor. 

0° 


Note 1. Sensors include multifrequency radiometers, boresight camera, 
analog recorder. 


2. Altitudes should be flown exactly as listed ±30 ft, all passes level 
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B. 


C. 


Recorded Flight Data 

1) Radiometer outputs and radiometer calibration data on PCM 
system and 4-pen analog recorder (use one sec integrators) 
(Record 4 frequencies plus time.) Be sure that in-flight 
calibrations are PCM recorded. 


2) 


3) 

4) 


5) 

6 ) 

7) 

8 ) 

9) 


Radiometer Frequencies 


Li -Band 
X-Band 
K-Band 
K^-Band 


1.42 GHz 
10.625 GHz 
22.2 GHz 
31.4 GHz 


Complete ADAS 

Boresight camera and RC8 imagery during all recording 
periods. (Must be time or frame correlated with radiometer 
data . ) 


Air temperature at each elevation. 

Air pressure at each elevation. 

Humidity at each elevation . 

Occurrence of water droplets (if any) at each elevation. 
Quantitative cloud cover - note height of cloud base and 
cloud tops m flight line (if any). 

Record laser altimeter data at applicable lower elevations. 


Ground Truth Data 


1) Nominal wind speed (steady or gustmg). 

2) Nominal wind direction. 

3) Nominal wave height. 

4) Nominal water temperature. 

5) Description of weather m general area. (Meteorological 
maps.) Presence of fronts, cloud cover descriptions, ceiling 
(if any). 
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6) Cold load and hot load calibrations should be run within a 
day or two "before" and "after" the flight, 

7) Types of materials traversed during passes {e.g., sand, 
grassland, trees, crops, etc.). 

The major difficulty in accomplishing the antenna pattern mission 
is m locating suitable targets. Ideally, a site should be seLected where 
there are a small number of well-defined terrain types. Within each type 
there must be negligible variation of the microwave properties of the sur- 
face, but from type to type there must be considerable variation Also, all 
but one of the terrain types must occur m patches whose sizes appear to be 
many tens of beamwidths wide at the lowest altitude at which the aircraft 
will fly, and from one to a few beamwidths wide at the highest altitude. The 
remaining terrain type would be the "background" within which the patches of 
the other terrain types lie. {It would probably be acceptable if there were no 
single background, as long as all of the patches met the size requirements 
stated above. ) 

The aircraft would then be flown over the site at high and low alti- 
tudes. The sidelobe correction procedure would be applied to the high alti- 
tude data, where sidelobe effects would be expected to introduce significant 
errors. The corrected results should then be the brightness temperatures 
of each of the terrain patches. To confirm this, the low altitude data would 
be used, because, mthe low altitude case, both the main lobe and the side- 
lobes would generally be seeing the same patch of terrain. Thus, except for 
"glitches" when crossing a boundary, there would be no significant sidelobe 
errors. 

It would, of course, be necessary to assure that there were no 
effects due solely to the diffei ent path lengths at the two altitudes (or, that 
such effects were adequately compensated for). This could be done by find- 
ing large terrain patches at the same site or at a nearby site. These 
patches would have to be many tens of beamwidths wide even at the high 
altitude. Flights over such patches would then produce data where errors 
due to sidelobes would be negligible from both high and low altitudes. Thus, 
any effects due to altitude itself would be isolated. 
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Antenna pattern flight test data were collected over various dikes 
m Galveston Bay, over a pond located m farmland area, and over the end of 
a cement runway. These data are less than ideal but can be used as a check 
on the validity of the antenna pattern correction scheme. 
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Section 5 


CONCLUSIONS 


Equations, flow diagrams and data reduction routines have 
been developed to correct MFMR data for atmospheric contributions and 
antenna pattern effects. The equations and techniques for performing these 
corrections are discussed m this report. Computer programs are being 
developed by Manned Spacecraft Center to apply these corrections to MFMR 
data. Equations and technical details for programming and debugging of 
the MFMR data correction routines were piovided as part of the data study. 

MFMR missions are flown for a variety of purposes and require- 
ments for atmospheric and antenna pattern corrections vary from mission 
to mission. Flexible data correction techniques have been developed to 
accommodate the range of experimental conditions The effectiveness 
of these correction techniques depends on proper mission planning and 
experimental procedures. MFMR mission planning must include 
consideration of data requirements and mission constraints for performing 
adequate corrections. 

Atmospheric contributions to MFMR measurements can be deter- 
mined by (1) computation of atmospheric attenuation and emission using 
suitable mathematical expressions and known or estimated meteorological 
conditions, and (2) direct MFMR measurements of atmospheric radiation 
at selected heights above the terrain. Each approach possesses advantages 
and limitations. Both are useful for MFMR data analysis. The first method 
consists of assuming that atmosphenc conditions are sufficiently well known 
(e. g , from data supplied by a nearby Weather Bureau station or by direct 
meteorological measurements from the aircraft) that an a priori calculation 
of the atmospheric effects can be made. The second approach assumes that 
the necessary parameters for an a priori calculation are not known (or 
known only partially) and thus requires that some measurements be per- 
formed before atmospheric corrections can be made. The latter approach 
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is more general, but has the concomitant disadvantage of requir- 
ing some flight time to be set aside specifically for atmospheric measure- 
ments. In addition, the calibration of the radiometers must be accurate 
for the second method to succeed. Both approaches are utilized m the 
microwave radiometer data study The reported equations and methods are 
sufficiently broad m scope for use with other microwave radiometers such as a 
C-Band system which is being considered for use m the Earth Resources 
Aircraft Program 

The effect of antenna patterns on MFMR measurements can be 
estimated from a knowledge of the antenna gam functions and the corres- 
poinding brightness temperature of each solid angle element where the 
antenna gam is significant The available patterns appear to be adequate 
for determining the total power m the sidelobes. Determination of the 
brightness temperatures for each significant solid angle element is made 
difficult by the ground scan which is along a single line. For most types 
of terrain, information will be needed on the terrain on all sides of the 
mam lobe. It is expected that temperature values accurate enough for this 
purpose can be estimated using the optical photography obtained at the same 
time as the microwave data. The antenna pattern correction can be time 
consuming and use of this correction is impractical for complex terrain 
data where a large catalog of brightness temperature signatures, is required 

Flight experiments and test procedures were developed and 
implemented to substantiate the equations and data reduction procedures 
for correcting MFMR data Experimental data taken to substantiate the 
atmospheric correction routine were inconsistent Unrealistic brightness 
temperatures were ottamed using antenna and radome loss values furnished 
for use m this study. Significant improvement was achieved by using an 
empirically determined loss factor. However, quantitative results could 
not be obtained due to the sensitivity of measurement data to the empirically 
determined loss factor. After an adequate calibration of the MFMR is 
achieved, verification tests of the atmospheric and antenna pattern routines 
can be completed The microwave radiometer data processing routines 
developed under Phase II of the data study will, for the first time, enable 

routine application of atmospheric and antenna pattern corrections. These 
refinements will facilitate use of MFMR data by the scientific community 
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